Stabilizing the Discrete Vortex of Topological Charge S = 2 
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We study the instability of the discrete vortex with topological charge S — 2 in a prototypical 
lattice model and observe its mediation through the central lattice site. Motivated by this finding, we 
analyze the model with the central site being inert. We identify analytically and observe numerically 
the existence of a range of linearly stable discrete vortices with S — 2 in the latter model. The range 
of stability is comparable to that of the recently observed experimentally S = 1 discrete vortex, 
suggesting the potential for observation of such higher charge discrete vortices. 



Introduction. In the past decade, lattice systems de- 
scribed by differential-difference equations in which the 
evolution variable is continuum and the spatial variables 
are discrete, have been a subject of increasing interest . 
These systems appear in many diverse physical contexts, 
describing, e.g., the spatial dynamics of optical beams 
in coupled waveguide arrays in nonlinear optics |2j, the 
temporal evolution of Bose-Einstein condensates (BECs) 
in optical lattices in soft-condensed matter physics [3|, 
the DNA double strand in biophysics [dj, and so on. 

One of the principal directions of interest in these lat- 
tice systems consists of the effort to understand the fea- 
tures of their localized, solitary wave solutions. In two 
dimensions, such structures can be discrete solitons 
or discrete vortices (i.e., structures that have topologi- 
cal charge over a discrete contour) |6(. In the past two 
years, there has been a considerable effort towards the 
observation of both entities in the context of optics, uti- 
lizing photorefractive crystals: regular discrete solitons, 
dipole solitons, soliton-trains, soliton-necklaces and vec- 
tor solitons were observed |7j , while two groups were inde- 
pendently able to experimentally produce robust discrete 
vortex states || . On the other hand, experimental devel- 
opments in the physics of BECs closely follow with promi- 
nent recent results, including the observation of bright, 
dark and gap solitons in quasi-one-dimensional settings 
, and with the generation of similar structures in higher 
dimensions appearing within experimental reach 10]. 

The above discussed experimental realization of dis- 
crete vortices of topological charge S = 1 || (i.e., with a 
2tt phase shift around a discrete contour) poses the ques- 
tion of whether higher topological charge discrete vortices 
could also potentially be experimentally realizable. The 
most natural higher topological charge state to consider 
is then the vortex with S = 2. However, lattice com- 
putations with a prototypical discrete model, namely the 
the discrete nonlinear Schrodinger (DNLS) equation, had 
identified that mode to be always unstable In this 
work we will revisit this topic and examine in some detail 
the instability of the S — 2 discrete vortex in the frame- 
work of the DNLS equation, which, in different variants, 
is relevant to all of the above mentioned, experimentally 
tractable settings. Our scope is to offer some insight on 



the nature of the instability, which will, in turn, allow us 
to suggest an explicit mechanism for its stabilization of 
this vortex by means of the inclusion of an impurity at its 
center. We analyze the latter case in detail and establish 
(analytically and numerically) the stability of the S = 2 
vortex in that setting for parametric regimes similar to 
the ones for which the discrete vortex of S — 1 has been 
found to be stable. 

Analytical Results. We consider the DNLS equation, 

i^n.m ^^2^n,m |^n,m| ^n,m; (1) 

where u is the complex field (the envelope of the elec- 
tric field in optics or the wavefunction in BECs), e 
is the coupling constant (the "tunneling rate" between 

Sites/wells), while A 2 ti„, m = W„+l,m+«n-l,m+«n,m+l + 
U n m — 1 

- A.u n , m is the discrete Laplacian. We seek 
standing wave, localized solutions in the standard form 
Un.m = 4 > n,m^ l ^~ A ^ t ■ Our approach is based on the 
Lyapunov-Schmidt theory for the existence of solutions 
[ll( and on linear stability analysis for tracing the stabil- 
ity eigenvalues of the corresponding solutions, similarly 
to what was done for discrete solitons and vortices in [12j . 

Our starting point is the so-called anti-continuum limit 
of e = where the nonlinear oscillators of our two- 
dimensional lattice are uncoupled. We excite a discrete 
vortex in that limit by choosing a contour containing 8 
sites ((-1,-1), (-1,0), (-1,1), (0,1), (1,1), (1,0), (1,-1), (0,- 
1)), where the corresponding solution is e %6n < m on these 
sites, while it is elsewhere (selecting without loss of 
generality // = 1). The motivation for the choice of 
the phases over this discrete contour is that we aim to 
construct a solution with S — 2 over the relevant con- 
tour, hence the real part of the configuration should 
behave as cos(2#), while the imaginary part as sin(20), 
which in turn immediately implies that we should choose 
0n,m — Jtt/2, where j = 1, . . . , 8 is an index over contour 
sites. Below, we briefly discuss the general theory that 
would apply to any solution over the relevant contour, 
and then focus on the discrete vortex with S = 2. The 
stationary state equation for 4> n ^ m is given by: 

1 2 \ i 



f(4>n,m, 4>n,m, e ) = (1 



e(A 2 +4)0„, m ,(2) 



and its complex conjugate f{4>n,m,4>n,m,t) = 0. The 
linearization operator for these two difference equations 
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,(3) 



with s n ' im >u nim = u n + n ',m+m'- Then, the solvability 
condition of the Lyapunov-Schmidt theory (allowing to 
continue a solution valid for £ = to e ^ 0) mandates 
that the projection of the eigenvectors of to the 

Eq. J5J and its conjugate is null. To 0(e), this con- 
dition provides the bifurcation function constraint gj = 
sm(6j - dj+i) + sin(9j - = 0, (for j = 1, . . . , 8 with 
periodic boundary conditions) that was algebraically ob- 
tained in |l4|. This is naturally satisfied for the S = 2 
vortex with 9j — 9j-\ — tt/2 discussed above. How- 
ever, computing the Jacobian matrix (A4i)j,fc = dgj /d6k 
of the bifurcation function g , one can observe that its 
eigenvalues are for the case of our S — 2 solution and, 
hence, second-order reductions are necessary to adjudi- 
cate on the existence/stability of the 5 = 2 vortex. 

Expanding = (f>° n m + €(j)} l m + 0(e 2 ), one can ob- 
tain the corresponding equations for the 0(e) correction 
to the solution profile 4> n as: 



(1 
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) 2 I 1 lim = (A 2 +4)< m . (4) 



The solution of Eq. (@J can be found as 
1 



n.m 



[cos( 



)+cos(0 i+ i -9j)]e^, (5) 
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discrete contour while 
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Using 1 to obtain the next order cor- 
the bifurcation function, we get: g 2 — 

'j+i - Qj) [ cos (8j - fy+i) + cos(9 j+2 - 9j+i)] + 

2 -Vj-l - ^)[ COS (%-^-l)+ COS (%-2-^-l)i + 

8ia(0j — 6j + 2)+sin(9j —^+4) (^,2+^,4 + ^3,6 +^,8); with 
1 < j < 8 and S denoting the Kronecker symbol. Once 
again the bifurcation condition is satisfied for our vortex 
of S = 2. However, the eigenvalues of the corresponding 
second-order Jacobian M2 are not identically zero 
and can be used to establish (in conjunction with the 
bifurcation condition being identically satisfied) the 
persistence of the vortex of topological charge S — 2 in 
the vicinity of e = 0. 

Furthermore, the Jacobian M2 of the second order re- 
ductions can be computed explicitly as: 



M 2 = 
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(6) 



By using the expansion 

' (7) 

one can study the stability of the discrete vortex of S = 2. 
Furthermore by expanding the eigenfunction in Taylor 
series in e (as we did above for the solution <fi) and cor- 
respondingly the eigenvalue A as A = eAi + 0(e 2 ), it can 
be shown (see ^1 f° r details) that the Jacobian M.2 can 
be directly connected with the eigenvalue correction Ai 
(for eigenvalues bifurcating from 0, which are the natu- 
ral sources of potential instability in the DNLS problem). 
The relevant equation connecting M.2 and Ai is the (re- 
duced, i.e., 8 x 8) eigenvalue problem of the form: 



M2C = Ai£ 2 c 



A? 



(8) 



where {C2)j,k — 1 if j = k — 1, — 1 if j = k + 1 and 
if \j — k\ 7^ 1. Using the discrete Fourier transform, one 
can obtain from JSJ the characteristic polynomial: 



A| - 2\i l-(-l) J -8sin 



2 



! sin 2 — x 
4 



(1 - (-l)^2sin 2 ^) 



0, 



j= 1,2,3,4, (9) 



which provides the leading order approximations to the 
eigenvalues of the S = 2 vortex as follows: in the 
neighborhood of A = 0, the vortex will have three 
eigenvalues of algebraic multiplicity four: A = and 
A = ±\/2ei, while it will have two simple imaginary 
eigenvalues A = ±\/ V80 + 8ei and two real eigenvalues 
A = ±\/ \/80 — 8e. Among the latter, the positive one 
is the reason for the S — 2 vortex being always (i.e., for 
any e 7^ 0) unstable, as was numerically observed 0]. 

The examination of the real eigenmode leading to the 
direct instability of the S — 2 vortex (that has support 
over the central, i.e., (0,0), site), as well as the appar- 
ent mediation (in numerical experiments-see below) of 
the instability by means of the central site, lead us to 
consider the possibility of having an "impurity" at the 
central site, e.g., a strong localized potential such as a 
laser beam in BECs or an inhomogeneity in the photore- 
fractive crystal, enforcing 0o,o = 0. In such a case, the 
bifurcation function g 2 lacks the last term (encompass- 
ing the Kronecker symbols), since these are interactions 
"mediated" by the now inert (0,0) site. Furthermore, the 
second order Jacobian is now much simpler and acquires 
the form: {M.2)j,k — 1 for j = k, — 1/2 for j = k ± 2 and 
for \j— k\ 7^ 0, 2. One can then repeat the calculation of 
the eigenvalues in the problem of Eq. (jHJ), via the discrete 
Fourier transform, to obtain the characteristic equation: 



Ai + 2i sin( — ) 
4 



= 0, 



3 = h . 



(10) 
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This results into three eigenvalues of algebraic multiplic- 
ity four, namely A = and A = ±eiv2/2. There are 
also two double eigenvalues A = ±2i. The crucial ob- 
servation, however, is that in this case, there are no real 
eigenvalues immediately present as e ^ and hence the 
discrete vortex with 5 = 2 will be linearly stable, due 
to the stabilizing effect of the impurity (or, to be more 
precise, due to the absence of the instability mediated by 
the (0,0) site). We now turn to numerical investigations 
to examine the validity of these findings. 

Numerical Results. We identify unit frequency solu- 
tions with topological charge 5 — 2, by initializing the 
exact solution configuration at the e = limit of Eq. 
(fT|) and then using continuation over e, combined with 
a contraction mapping for the solution of the nonlinear 
system of (algebraic) equations to identify the exact (up 
to a prescribed numerical accuracy) numerical discrete 
vortex. We then perform linear stability analysis, using 
the expansion of Eq. (0, to obtain the eigenvalues A, 
and their corresponding eigenvectors. 

Figure ^ shows a typical example of the discrete vor- 
tex for e = 0.2 in the regular DNLS model. The middle 
and right panels show the real and imaginary part of the 
solution, clearly emulating cos(26*) and sin(2#) over the 
lattice contour of interest. The linear stability analysis 
of this vortex is shown in Fig. [21 One can observe that 
both for the imaginary (top left panel), as well as for 
the real (bottom left panel) eigenvalues, the predictions 
of the perturbation theory (dashed line) are extremely 
accurate in comparison with the full numerical results 
even for e up to 0.25. Clearly as e increases, higher order 
phenomena become relevant such as the splitting of the 
quartet of eigenvalues at A = ±^/2ei, or the collision of 
the simple pair of eigenvalues with A = ±y/ \/80 + 8ei 
with the bottom edge of the continuous spectrum (which 
is at A = ±i), resulting into a Hamiltonian Hopf bifurca- 
tion and a complex quartet of eigenvalues for e > 0.23. 
Additional such quartets appear for larger values of e. 
However, the solution is always unstable due the real 
eigenvalue pair A = ±y \/80 — 8e, whose eigenfunction is 
shown in the right panel of the figure. Notice that the 
latter has support over the central site, predisposing us 
for the role of this site in the instability development. 

The corresponding predictions/numerical results for 
the model with the impurity (i.e., with (0,0) inert) are 
shown in Fig. [3J The top panel illustrates the eigenvalue 
of multiplicity four with A = ±ei\/2/2 and with multi- 
plicity two A = ±2i, which are again in excellent agree- 
ment with the numerical findings. The real part clearly 
indicates the absence of an instability for small e. Such 
an instability arises due to collision of the eigenvalue pair 
with the continuous spectrum and is present for e > 0.36. 
It is crucial to note here that the 5 = 1 vortex was found 
to be stable for e < 0.38 in This illustrates that 

the present mechanism stabilizes the 5 = 2 vortex for a 




FIG. 1: The discrete vortex is shown for e = 0.2. The left 
panel shows contour plots of the real (top) and imaginary 
(bottom) parts. The middle panel shows a 3d rendering of 
the real part, while the right panel shows a similar 3d plot of 
the imaginary part of the vortex. 




FIG. 2: Linear stability analysis of the S = 2 vortex. The left 
panel shows the imaginary (top) and the real (bottom) part 
of the point spectrum eigenvalues bifurcating from A = 0. 
The solid line shows the numerical result, while the dashed 
line indicates the theoretical prediction (see text). The most 
unstable eigenmode pertaining to the real eigenvalue for e = 
0.2 is shown (real part: top; imaginary part: bottom) in the 
right panel. 

parametric region comparable to that of the 5 = 1 vor- 
tex, hinting that it could be experimentally feasible to 
trace such a configuration similarly to what was done for 
the 5 = 1 case |8j. The right panel shows the real and 
imaginary part of an unstable eigenmode for e = 0.4. 

Finally, to examine the dynamical development of the 
instability and to compare/contrast the dynamical fea- 
tures of the two models (in the absence and presence, 
respectively, of the impurity), we have conducted direct 
numerical experiments. The main results are shown in 
Fig. 0| for two representative cases (namely e = 0.2, 
where the former case is unstable, while the latter is sta- 
ble, and e = 0.4, where both models have unstable 5 = 2 
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FIG. 3: Same as in Fig. 2, but for the model with the impurity 
(i.e., with (0,0) inert). The instability is absent in this case. 
The right panel shows the principal eigenmode of instability 
for e = 0.4. 
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FIG. 4: For identical initial conditions (see text), e = 0.2 
(left) and e = 0.4 (right), the solution is shown in the relevant 
central sites as a function of time, in the absence (top) and 
presence (bottom) of the impurity. The central site (dash- 
dotted line) and the sites (— 1, 0) (solid line) and (1, 0) (dashed 
line) are shown (the latter two are the only ones that remain 
excited in the configuration for long times). 



ity interval of the experimentally observable S = 1 state. 
Numerical findings conclusively corroborate this picture 
both at the level of linear stability analysis (found to be 
in excellent agreement with the theoretical predictions) 
and at the one of direct numerical experiments. We be- 
lieve that this opens the path for observation of higher 
charge discrete vortices and renders experimental work 
in this direction particularly timely. 

Numerous constructive discussions and comments by 
D.E. Pelinovsky and partial support by NSF-DMS- 
0204585, NSF-CAREER, and the Eppley Foundation for 
Research (PGK) are gratefully acknowledged. 



vortices). In both cases, we have simulated both models 
up to t — 200, initializing them with identical initial con- 
ditions consisting of the vortex with a perturbation (mul- 
tiplied by 10~ 4 ) in the eigendirection of the right panel 
of Fig. 0for e = 0.2 and of Fig. [3 for e = 0.4. In the case 
of e = 0.2, we observe that the DNLS vortex becomes 
unstable, whereas in the presence of the impurity the so- 
lution is completely stable (exhibiting oscillations at the 
order of the initial perturbation). The instability for the 
DNLS vortex appears to be mediated by the central site 
(dash-dotted) line, which eventually settles at a rather 
small amplitude. The final configuration finds 6 of the 
8 (initially) participating sites at the vortex with near- 
zero amplitudes, while only two sites (shown by solid 
and dashed line) remain excited in an asymmetric con- 
figuration with a long-lived breathing (weak) exchange of 
power between them, mediated principally by the central 
site. The instability sets in for t rs 45 in this case. For 
e = 0.4, the regular DNLS becomes unstable even faster 
(f ps 30) and the dynamics is similar. For the case with 
the impurity the instability sets in at much longer times 
(t ps 80), as expected by the much smaller value of the 
corresponding principal eigenvalue real part. Further- 
more, while only two sites remain excited in this case as 
well, the configuration is no longer asymmetric and the 
oscillation of power can be identified (data not shown) as 
being caused by the small amplitude exchange of power 
around the vortex contour (recall that in this case the 
central site is in this case inert). 

Conclusions. In this work we have revisited the topic 
of discrete vortices of topological charge 5 = 2. We 
have explicitly discussed and illustrated their instability 
in the prototypical lattice model of the discrete nonlinear 
Schrodinger equation and have traced its source in the ex- 
change of power made available through the central site 
of the vortex. We have thus proposed to consider a model 
with an impurity (an inert) site at the center of the vor- 
tex. Examination of the stability problem in the latter 
case shows the absence of linear instability for a regime of 
coupling strengths comparable to that of the linear stabil- 
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